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Abstract 



We examine a new application of the Holstein-Primakoff realization of the simple har- 
monic oscillator Hamiltonian. This involves the use of infinite-dimensional representations 
of the Lie algebra su{2). The representations contain nonstandard raising and lowering 
operators, which are nonlinearly related to the standard and a. The new operators 
also give rise to a natural family of two-oscillator couplings. These nonlinear couplings 
are not generally self-adjoint, but their low-energy limits are self-adjoint, exactly solvable, 
and stable. We discuss the structure of a theory involving these couplings. Such a theory 
might have as its ultra-low-energy limit a Lorentz- violating Abelian gauge theory, and we 
discuss the extremely strong astrophysical constraints on such a model. 
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1 Introduction 



The simple harmonic oscillator is one of the best understood systems in quantum mechan- 
ics, with applications in essentially all areas of physics. However, there may remain many 
interesting properties of this system that have not been fully understood or elucidated 
(see PP, for instance). In this paper, we present an example of this. This example provides 
a possible way for nonlocal interactions to stabilize a Lorentz-violating modification to 
the free photon sector of the standard model. 

We shall show how the Holstein-Primakoff realization of the angular momentum Lie 
algebra su{2) may be used in connection with the simple harmonic oscillator We first 
review the free oscillator, considered using the infinite-dimensional representations of this 
algebra. This leads naturally to a study of new nonlinear couplings between multiple 
identical harmonic oscillators. These couplings can be viewed in one of two ways. Taken 
at face value, the interactions are not self-adjoint, and the energy eigenvalues need not 
be real; however, the Hamiltonians in question do possess low-energy limits which are 
self-adjoint. The second possible viewpoint would be to take the formula which defines 
the real eigenvalues in the low-energy regime and extend that formula to cover the full 
range of the quantum numbers. This has the advantage of ensuring unitary evolution, 
but the additional complexity of the Hamiltonian is a corresponding disadvantage. 

We shall describe how these new interactions may manifest themselves in a physical 
theory, with particular emphasis on couphngs between the two polarization modes of the 
free electromagnetic field. However they are interpreted, the interactions we shall discuss 
may be relevant as part of a nonlocal, Lorentz-violating quantum field theory. Recently, 
there has been a great deal of interest in the possibility that Lorentz symmetry may not be 
exact in nature. A violation of this kind of fundamental symmetry could arise as part of 
the novel physics of the Planck scale. Relics of this violation would then persist even in the 
low-energy effective theory. The general local Lorentz-violating standard model extension 
(SME) has been developed and the stability and renormalizability [7j of this 

extension have been studied. Lorentz violation is a very interesting area of theoretical 
physics, because even superficially simple questions about its physics may have subtle 
and even ambiguous answers. For example, the study of the gauge invariance properties 
of and finite radiative corrections to Lorentz-violating field theories has proven to be a 
fruitful source for new theoretical insights jHl El UDl UH 1121 UH]- 

One significant difficulty with Lorentz-violating quantum field theories is that they 
frequently exhibit problems with stability [UlEj- Yet it has been suggested that some form 
of nonlocality might overcome this problem jHlE]. The harmonic oscillator interactions 
we consider can provide a concrete example for how this stabilization might work, if we 
adopt the second interpretation of these interactions as described above. In the limit 
of very low energies, the nonlocal interactions may couple together the two polarization 
modes of a free photon in exactly the same way as would a local, renormalizable operator 
from the SME P[ Ej • However, the nonlocality ensures that the theory remains stable. 
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even for very large photon numbers. Weak forms on nonlocality have also been considered 
in other Lorentz- violating contexts |16j . 

The kinds of infinite-dimensional representations of su{2) that we will consider have 
also been introduced in the context of the Dirac Coulomb problem [T71 ITH] and in gen- 
eralizations of the Dirac monopole where they provide useful insights. Moreover, 
other Lie algebra also possess infinite-dimensional representations that are not represen- 
tations of any corresponding Lie group. These representations might be useful in the 
study of certain quantum-mechanical systems, through a generalization of the techniques 
used in jTTl UHl HH] or developed in this paper. 

2 Single-Oscillator Operators 

We shall begin by reviewing the Holstein-Primakoff realization of the harmonic oscillator 
raising and lowering algebra. Most frequently, when one studies the representations of 
su{2) in connection with quantum mechanics, one is interested only in the finite-dimen- 
sional representations, which are countable and parameterized by the total angular mo- 
mentum i. One considers an operator J = (Ji, J2, J3) = (|^+ + — ^J-i -^s); 
with standard commutation relations. Beginning from a highest weight state with 
J+\i) = 0, one constructs each of the 2i+l states \i — s) by acting on \i) with J_ s times 
and normalizing appropriately. 

When 2i is a nonnegative integer, there are only these 2i + l states, because J_| — ^) = 
0. However, for more general values of the highest weight, the sequence of states does not 
terminate. Instead, one constructs an infinite tower of equally spaced states. This tower 
of states has a structure which is identical to that of the simple harmonic oscillator. 

In fact, if we begin with a highest weight state |A), with J3IA) = A|A) (where 2A is 
not a nonnegative integer), we may construct states \X — n) for all nonnegative integers 
n, using J_|A — = \J X{X + 1) — (A — n){X — n — 1)|A — n — 1). Then the Hamiltonian 

Hx = —ioJ^ + (^A + has nondegenerate eigenvalues (n + lo, which are precisely 
the energy eigenvalues of a harmonic oscillator (when we set h = 1). Since a quantum- 
mechanical system is entirely specified by its Hilbert space and the Hamiltonian acting on 
that space, this is equivalent to an alternate description of the harmonic oscillator. (Note 
that this description is completely distinct from Schwinger's development of the angular 
momentum algebra in terms of harmonic oscillator states [201 ■) 

We should point out that the state space is a representation only of the Lie algebra 
su{2), not of the Lie group SU{2). That this is the case should be clear from an ex- 
amination of the spectrum of the Hamiltonian Hf^ = ujfi ■ J. This would be a Zeeman 
effect Hamiltonian if the representation of J were finite-dimensional, with J transforming 
in the adjoint representation of SU(2). The SU{2) symmetry would then dictate that 
the eigenenergies should be independent of the direction of n. However, those energies 
are clearly not independent of n for the infinite-dimensional representations we are now 
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considering. Specifically, if h = —63, then we have a harmonic oscillator system, but if 
n = +63, then the energy is not bounded from below. Therefore, the infinite-dimensional 
operators cannot form a representation of the group SU{2). 

We shall now re-label our states, so that they match the usual harmonic oscillator 
nomenclature. We make the replacement |A — n) — *■ |?7,). Now, in addition to the standard 
harmonic oscillator raising and lowering operators and a, we have a new set of raising 
and lowering operators J_ and J+, given by 



The square root operators are to be interpreted as having eigenvalues equal to the square 
roots of the eigenvalues of the operators inside, with the same eigenvectors. For a me- 
chanical oscillator, we may further express the J± in terms of the position and momentum 
operators x and p, using the usual linear relations connecting x and p to a'' and a; however, 
this substitution must be made with the understanding that the proper interpretation of 
the J± operators requires the use of the discrete "number of quanta present" basis of 
states. 

The J± give us a new family of raising and lowering operators, parameterized by A. 
These operators are distinct from the and a for all finite A. This is clear from the 
differing commutation relations [a,a^] = 1 and [J+, J_] = 2J3. In general, for large n 
(large in comparison with A and unity), the matrix elements of J± are larger that those 
of a and by a factor of 0{^/n). However, the new operators do include the a and as 
limiting cases. As A ^ 00, —>■ a and —>■ . (These limits are to be interpreted 
in terms of the matrix elements of the operators involved, and all half-integral values of 
A must be avoided as the limit is taken.) 

We can see from ^ and that J+ 7^ j1 for finite, non-half-integral A, because the 
square roots in (jH)) and may become imaginary. This means that Ji and J2 are not 
self-adjoint, a difficulty which we glossed over when we discussed the Hamiltonian Hf^. 
The Casimir operator is self-adjoint, however. Moreover, in the basis of eigenstates 
of J3, J_|_ = J^. Since the matrix elements of J_ and differ only by phase factors 
in the J3 basis, there is no ambiguity in defining the entire state space starting from the 
ground state. Finally, we point out that if A is large and positive, then the non-self-adjoint 
character of J does not become apparent unless n is at least comparable to 2 A. These 
facts will prove important when we discuss the coupling between two harmonic oscillators. 



J_|n) 
J+|n) 




(1) 
(2) 



or 




(3) 
(4) 
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3 Multiple Coupled Oscillators 



We shall now consider a novel application of this description of the harmonic oscillator. 
We may determine the energy eigenvalues exactly for certain systems of coupled identical 
oscillators in which the couplings are nonlinear. Let us consider the Hamiltonian 

H = H^ + H^ + gH,,,, (5) 



H 



-uJaz + ( a + i ) ct; 



+ 



-ujJbz + (^/i + w 



+ gHir^f (6) 



J A and Jb are two independent vectors of operators of the type we have been considering, 
corresponding to the highest weights A and yU, respectively, i^int is an interaction, whose 
form we shall discuss shortly. This Hamiltonian has three adjustable parameters. A and 
IX determine the structure of the harmonic oscillator representations that we are using. 
However, like g, they may be seen as parameters describing the interaction, because we 
have shifted the total energy in such a way as to make the spectrum of the free oscillators' 
Hamiltonian Hx + H^ independent of both A and fi. These two parameters may be chosen 
freely, subject to the condition that neither 2A nor 2/i is a nonnegative integer. 

We shall choose an interaction gH^^i that is similar in form to gJ^ ■ Jb- If 0, 

— * — * 

gJA • Jb is not self-adjoint; however, if A and // are large and positive, this problem will 
not be apparent in the vicinity of the ground state. So we take iJint to agree with Ja ■ Jb 
when Ha + ns + 1 < min(2A, 2fi), where ua and ub are the principal quantum numbers 
of the two oscillators. We may then consider the effects of this interaction within this 
restricted ("low-energy") regime. 

We may solve the restricted Hamiltonian using the ordinary techniques for the addi- 

^ — * — * — * 

tion of angular momenta. The operator J = Ja + Jb has highest weights of the form 
\ + fi — i for all nonnegative integers i, and each value of i corresponds to a single ir- 
reducible component of the representation. The Clebsch-Gordon coefficients for these 
representations can be calculated by the standard method of applying lowering operators 
and using Gram- Schmidt orthonormalization. However, our primary interest is in the 
energy levels. 

The eigenvalues of Hx + are just {jia + ^b + 1)ciJ. When we change the basis, to 
use the "total angular momentum" J, this part of the Hamiltonian becomes 

Hx + Hf,^ -a; Jg + (A + + l)oj. (7) 

The eigenstates of the system are parameterized by the highest weight \ + fi — i and by 
the "number of quanta present" (that is, the number of applications of J_ on the highest 
weight state required to produce a given state), n. If we denote these states by \i, n), then 
it is clear that {Hx + Hn)\i, n) — {i-\-n-\- l)a;|i, n). Since i and n have the same range as 
Ha and tlb (all must be nonnegative integers), this verifies that the free system has the 
same spectrum in each basis. 
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The coupling term must be calculated in the J basis. This is easily done, using 
Ja-Jb = \ {J^ ^ Ja ~ Jb)- For a state of specified i, Ja ■ Jb has the eigenvalue Xjj — 

(X + H+^)i + \i'^. So the total energy is 



Ei^n = uj{i + n + 1) + g 



A/i - ( A + + ^) « + V 



This formula holds exactly in the entire low-energy subspace. 

We must now turn our attention to the general definition of ifmt- There are two 
natural ways to define this interaction. The first possibility is that ifjnt = Ja- Jb exactly, 
and the fundamental Hamiltonian is not self-adjoint. The second possibility involves a less 
drastic modification of the structure of the theory. We simply take the exact Hamiltonian 
to be defined by its eigenvalues, which have the form (jH)). There is a complete set of 
states {\i, n)} corresponding to these eigenenergies, and for i + n+ 1 < min(2A, 2/i), these 
states agree with the ones found above. The new Hilbert space is not equivalent to the 
old, two-oscillator Hilbert space, but the restricted, low-energy subspaces of these two 
Hilbert spaces are isomorphic. We shall henceforth adopt this second definition of ifmt 
(although most of our statements will concern only the low-energy subspace, where the 
two definitions are equivalent). 

A few words about the various parameters are now in order. Our coupled system is 
stable only if g > 0; if g < 0, the energy is not bounded below. We also see that the 
structure of the energy levels depends only on g and A -|- yu; the separate values of A and 
fi only affect the zero-point energy. 

This is a strongly nonlinearly coupled system. We may recover the more usual result 
for the spectrum in the presence of a linear coupling between the oscillators by setting 
X = fi and expressing and a in terms of the A — oo limits of J+ and J_. (Note that 
as A — s> oo, there are no problems with the operators being self-adjoint.) The interaction 
g (aAdB + '^a'^'b) is the limit of j (^Ja ■ Jb — Ja^Jb^ as A approaches infinity. The two 

terms ^Ja'Jb and ^JasJbs each commute with the noninteracting Hamiltonian; and, 
although they do not commute with each other for finite values of A, they do commute in 
the infinite limit. We may see this by evaluating the two operators for finite A in different 
bases. In the J basis, the first term is diagonal, with eigenvalues g X — (2 + i + , 
just as calculated above. The second term is diagonal in the Ja-, Jb basis, with eigenvalues 
|(A — nA){X — ub). As a ^ oo, these eigenvalues become g{X — 2i) and g{X — ua — ub), 
respectively. However, as we saw when we discussed the noninteracting case, the total 
number operator UA+ns = n+i is diagonal in both bases, so (^(A— n^ — n^) = g[X — n — i). 
Then the total energy shift, which is now the difference between g{X — 2i) and g{X — n — i), 
becomes simply g{n — i). That is, we have two decoupled oscillators with frequencies uj±g, 
which is just the usual result. 

More general interactions are also possible. We may replace gJA ■ Jb with an arbitrary 
function of Ja ■ Jb- Any such interaction will still commute with H\ + H^, and all the 
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same considerations will still apply. Other generalizations are possible as well, through 
the use of other well-known properties of su{2). For example, we may generalize to the 
coupling of identical oscillators, with the identity 



l<C<D<Af 




(9) 



However, it is important to keep in mind that many approximations that are typically 
used when one studies more than two interacting angular momenta will break down when 
working with infinite-dimensional representations of su{2). 



One of the simplest situations in which pairs of identical harmonic oscillators arise is 
in an Abelian gauge theory in 3-1-1 dimensions. An interaction of the form (jH)) might be 
relevant as a modification of the photon sector of quantum electrodynamics. For the gauge 
sector alone, this interaction could be introduced separately at each value of the photon 
momentum. The selection of a specific basis of polarization states for each momentum, 
and the assignment of the couplings A, fi, and g generally breaks Lorentz symmetry and 
may also break parity invariance. Moreover, the operators J± are nonlocal, since they 
involve the total energy present in a given mode of the electromagnetic field, and it 
is not possible to express this sort of interaction conveniently in terms of the ordinary 
electromagnetic field operators A'^ and F^'^. However, if we are willing to allow these 
modifications to the structure of the theory, the interaction is (in the absence of charges) 
exactly solvable. Since there is current interest in exotic modifications of QED, this type 
of interaction may be worthy of further investigation. 

We must also say a word about the gauge invariance of this QED modification. The 
Ja- Jb interaction is formulated in terms of creation and annihilation operators (i.e. in the 
canonical quantization formalism) , and so any discussion of gauge invariance will necessar- 
ily be complicated by the difficulties that are associated with the canonical quantization 
of gauge fields. However, at low energies, the interactions we have introduced are clearly 
consistent with gauge invariance in the following sense. If we specialize to the Coulomb 
gauge and quantize the transverse modes of the theory, then the low-energy interaction 
may be introduced without difficulty. It does not affect the number of polarization states 
(and is, in fact, strongly dependent upon this number), so it does not spoil gauge in- 
variance in this fashion. However, it is possible that it may damage gauge invariance at 
higher energies or when interactions with matter are considered. 

Although the interactions we are considering cannot be expressed simply in terms of 
the standard electromagnetic field operators, the ultra-low-energy limit of our theory could 
well be expressible in such a form. We shall shortly show that this is indeed the case for a 
particular class of models. The embedding of a ultra-low-energy Lorentz- violating effective 
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field theory within the framework of our fundamental theory is attractive for several 
reasons. First, the interactions we have considered have eigenenergies that are exactly 
known. Second, although Lorentz-violating theories may exhibit stability problems, our 
theory does not; the term in (jH)) ensures this. This example shows that there can 
exist stable nonlocal interactions which have local, Lorentz-violating Lagrangian theories 
as their ultra-low-energy limits. Third, such an embedding demonstrates the intriguing 
possibility that our conventional basis of polarization states many be inadequate for the 
description of the full Hilbert space of a more fundamental theory. 

We shall therefore examine the ultra-low-energy behavior of our QED modification, to 
see to what sort of effective theory it may correspond. Since the energy-level differences 
depend only on A + yU, we shall set A = /i. We take A to be a very large number, large 
in comparison with any relevant photon occupation number. (This is what we mean 
by "ultra- low-energy.") We may then neglect a^a compared to A in and It 
immediately follows that the energy eigenstates are approximately given by 

\i,n) ^ ^ ^ ^ ^ |0,0). (10) 



(^!)(n!) 

Deviations from this expression are suppressed by a factor of (9(A^^/^). The corresponding 
energies, in the same approximation, are 

Ei^n^ {n + i)uj + 2\gi, (11) 

where we have dropped the zero-point energy. 

These results imply that the polarization modes corresponding to = ^ {a\ + a)^ 

and at = {a\ — a^^ have different frequencies. The two frequencies are shifted from 
their mean value by ±^ ~ ±fi'A, and the effective Hamiltonian is 



ifeflf = (cu + 2Xg)a}_a^ + uoa^aj^. (12) 

This same effective Hamiltonian arises naturally in the context of a CPT-even, Lorentz- 
violating modification of the photon sector. For a theory with Lagrange density 



C = - -^Kr^.F^^F'^'' (13) 



(where fcrer/^i/ has the symmetries of a Riemann tensor and is double traceless), the ex- 
pressions for the photon modes' frequencies are (to leading order in ki^rf^u) (HI Ul IHj 

u;i = {l+p±a)\p\. (14) 

Here, p is the photons' 3-momentum, p = — and = ^ka/3k°'^ — p^, with ka/3 = 
ka^pyP^p" and p^ = {l,p/\p\). The approximate frequencies given by (fT^ and (fT^ cor- 
respond if p = a = gX. So any theory with p = a > for all p will reproduce the entire 
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low-energy behavior of our modified theory. Theories with this property indeed exist; for 
example, if kf^rfiu has the form 




(15) 




then p = a = —w'^, where = v^{u ■ p) — u^{v ■ p). If f = (V, Oj and = (0, u), with 

\u\ = 1, then p = sin^ 9, where 9 is the angle between u and p. This particular model 
is parity-preserving, with three independent parameters. 

Moreover, any theory with nonvanishing a will demonstrate a splitting between po- 
larization modes, as would arise in our s'u(2)-modified QED. There is therefore a large 
theoretical parameter space in which the theory given by ()13|) can be embedded in a su{2) 
coupling model. 

These embeddings all require that, for a fixed direction p, the coupling gX must be 
proportional to \p\. It would seem most natural for A, which represents the number of 
photons that must be present in a mode of the electromagnetic field in order for the failure 
of the polarization state basis to be apparent, to remain large for all values of \p\. We 
therefore speculate that A may be a |p|-independent (or ^-independent) constant, while g 
scales with the magnitude of p. 

Based upon astrophysical experiments, the physical value of a is strongly constrained, 
to parts in 10^^ or better [^|22]- This represents an even stronger constraint on the su{2) 
model, because A is necessarily very large. It follows that g = a/\ is correspondingly 
smaller. Our modification of QED is thus physically reasonable only in a very, very small 
region of parameter space. 

5 Conclusions 

In summary, we have studied an applicaiotn of the Holstein-Primakoff realization of the 
simple harmonic oscillator operator algebra. This has involved the introduction of a set of 
raising and lowering operators that obey the angular momentum commutation relations. 
The properties of these operators allow us to solve exactly for the low-energy behavior of 
a theory with a particular nonlinear interaction. The high-energy extension of this model 
may involve either a non- self-adjoint Hamiltonian or a new basis of states. 

This model has an interesting application in Lorentz-violating physics, where stability 
is typically a problem. The ultra-low-energy limit of the interaction we have considered 
resembles the effect of a small Lorentz-violating correction to free QED. At higher energies, 
the nonlocality of the special harmonic oscillator interactions we are considering can 
stabilize the Lorentz-violating theory. Physically, these kinds of effects are extremely 
strongly constrained by astronomical observations. However, this theory provides an 
useful insight into how nonlocality and Lorentz violation may combine to form a well- 
behaved theory. 
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